Spin-orbit coupling effects on quantum transport in lateral semiconductor dots 
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The effects of interplay between spin-orbit coupling and Zeeman splitting on weak localisation and 
universal conductance fluctuations in lateral semiconductor quantum dots are analysed: All possible 
symmetry classes of corresponding random matrix theories are listed and crossovers between them 
achievable by sweeping magnetic field and changing the dot parameters are described. We also 
suggest experiments to measure the spin-orbit coupling constants. 
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The effects of spin-orbit (SO) coupling on transport 
phenomena in chaotic quantum dots |l|,|2j recently at- 
tracted attention. Motivated by a puzzling modification 
of the variance of the mesoscopic conductance fluctua- 
tions with applied in-plane magnetic field, Halperin et. 
al. M suggested that the specific form of the spin-orbit 
interaction in a 2D electron gas based on semiconduc- 
tor heterostructures may be responsible for a series of 
crossovers not considered in the existing literature pi. 
It has also been noticed Hpf that spin relaxation in a 
quantum dot may be facilitated by the Zeeman field. 

The goal of this paper is two-fold: (i) we provide a 
complete quantitative theory for the interplay between 
spin-orbit coupling effects and Zeeman splitting in zero- 
dimensional semiconducting systems, identify all the pos- 
sible symmetries, and describe all the physically achiev- 
able crossovers; and (ii) we show that the SO cou- 
pling effects depend on the magnetic field orientation 
in anisotropic dots and discuss possible experiments en- 
abling one to measure directly the ratio between two in- 
dependent SO constants. 

The single-particle Hamiltonian of the system, H = 
Hq + u(r) is the sum of the free-electron dispersion term 
and a potential, u(r), consisting of a confining potential 
and a random potential of impurities. The free-electron 
term includes spin-orbit coupling, as a combination of a 
Rashba term and a crystalline anisotropy term (specified 
for (001) plane of GaAs), and Zeeman splitting energy 
due to the in-plane magnetic field @], B = IB, 
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where p = P — eA is the kinetic momentum, with P being 
the canonical momentum and A being the vector poten- 
tial describing the orbital effect of the magnetic field. 
Since (001) plane of GaAs has the symmetry of a square 
without inversion centre, C<i v , we choose the coordinate 
system (3:1,3:2) with axes along crystallographic direc- 
tions e*i = [110] and e*2 = [110] and rewrite Hq as 
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where Aj~ 2 = ot±g characterise the length scale associated 
with the strength of the spin-orbit coupling for electrons 
moving along principal crystallographic directions [ci,2,3 
are Pauli matrices, cr 2 — — crj, o\$ — of 3 ]. 

The Hamiltonian H = Ho + u(r) describes the electron 
motion in a lateral semiconductor dot coupled to metal- 
lic leads via two contacts, / and r, each with Ni, r > 1 
open orbital channels. Below we focus on the universal 
0D description applicable if B 



7, e z < E T ; L12 < A12. 
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Here, 7 = (Ni + N r )A/2irh stands for the escape rate 
into the leads, A = 2irh 2 /mA is the mean level spacing 
in a dot with area A ~ L1L2, and Et is the conventional 
Thouless energy. 

Our purpose now is to identify Hamiltonian (m) with an 
appropriate random matrix ensemble. Doing it directly, 
however, is not convenient. The reason for this is that 
on shell matrix elements of the velocity vanish due to the 
gauge invariance (the importance of this fact for SO in- 
teraction in quantum dots was first noticed in Ref. 0). It 
means that if the spin remained fixed during the motion 
of the electron, the effect of SO coupling would be just a 
homogeneous shift of the momentum space which could 
not change observables. To get rid of such terms fixed by 
gauge invariance, we perform the unitary transformation 
of the Hamiltonian as H — > H = U^HU with 



U — exp 
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Using the condition £1,2/^1,2 < 1, we expand H up 
to the second order in the coordinates and obtain 



H 



1 (-> A - CTz 

p — eA — a 1 — 
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A= B z [fxn z ]/2c; a± = [rxn z ]/2AiA 2 ; 



Y + h {0) + h {l) + u(r) 
(4) 
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Equation (^) indicates that the effects of SO coupling in 
the leading order at ez = 0, and of the orbital magnetic 
field are somewhat similar. This similarity is not a coin- 
cidence — in the leading order, the direction of the spin 
follows the motion of the electron: for an electron moving 
along a closed path, its spin spans the closed path too. 
Due to the motion in spin space, an electron accumulates 
extra Berry phase equal to the solid angle spanned by the 
spin. At weak SO coupling, this area is proportional to 
the geometrical area encircled by the electron path in 
the coordinate space resulting in an effect similar to that 
of Aharonov-Bohm flux. This analogy may be put on a 
quantitative level by noticing that the two energy scales 
characterizing both effects 
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have the same dependence on the shape and the disor- 
der in the sample. Here, k is the coefficient dependent 
on the geometry and A is the area of the dot. Random 
quantities M a p are the non-diagonal matrix elements of 
the magnetic moment of the electron in the dot. 



Term (0) is higher order in the SO coupling constant. 
However, it has a different symmetry from a±, therefore, 
its retention is legitimate. Its physical significance is to 
provide the spin-flips and, thus, the complete spin re- 
laxation, in contrast to a± which preserves correlations 
between spin up and spin down states. Quantitatively, 
the effect of 3\\ is characterized by the scale 
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The effect of the in-plane magnetic field is described 
by Eq. m). It includes the homogenous Zeeman split- 
ting /i^and the combined effect of the SO interaction 
and Zeeman splitting described by h^\ The latter can 
be envisaged as a deflection of the effective magnetic field 
from the direction given by external B, and it results in 
the spin relaxation associated with the energy scale 
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where x 12 are the non-diagonal matrix elements of the 
dipole moment of the electron in the dot. Quantity Sy 
depends on the geometry and on the disorder in the dot 
and may be estimated as S ~ AL 2 /Et, so that t\ <C ez- 
A similar energy scale has appeared in recent publica- 
tions P|0]j however, the symmetry of the corresponding 
term hy> was not indentified. 



TABLE I. Symmetries of the system in the absense of orbital magnetic field effect, tb 3> T e , 
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Having derived a Hamiltonian free of the gauge invari- 
ance constraints on the values of its matrix elements, we 
identify the symmetries of all relavant limits. The results 
are summarized in Tables 1 and O depending on the or- 
bital effect of the magnetic field. In these tables, the 
conventional parameter j3 describes time-reversal sym- 
metry of the orbital motion, s is the Kramers degeneracy 



parameter, and E is an additional parameter character- 
ising the mixing of states with different spins for strong 
Zeeman splitting. Parameters /3, E, and s completely 
characterize the statistical properties of the transport 
through the system as well as spectral correlations of the 
isolated dot. The straightforward generalization of the 
known results, 0], gives the following description of the 



two-terminal conductance of the dot (measured in units 
of ^jr) connected to the leads by reflectionless contacts 
with Ni and N r orbital channels: 
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The parametric dependences of the transport coeffi- 
cients in Eqs. (ll|lj) can be envisaged as a sequence of 
crossovers shown below as a function of Zeeman splitting 
energy and escape rate (increasing along the horizontal 
and vertical axes, respectively): 



l[lu] -► 4 [4m] -> 7[6u] 

T T T 7 

2[2u] -► 5[5u] -> 8[6u] fr 

T T / => ez 

3[3u] -> 6 [6m] 
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For a large number of channels, 1 <C Ni + N r <C Et/A, 
the result of Eqs. (0,03) can be simplihed: 
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and the crossover in WL behavior between various sym- 
metry classes in the Tables |£l| sketched in (O) can be 
described quantitatively using a diagrammatic perturba- 
tion theory J7|, as we outline below for the weak locali- 
sation (WL) correction. 



TABLE II. Symmetries of the system in the presense of orbital magnetic field effect, tb <C T e , 
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For a spin- 2 particle in a quantum dot with ballis- 
tic adiabatic contacts, the WL correction to the conduc- 
tance can be related to the lowest-lying modes of Cooper- 
ons in a singlet (L = 0), P<S° = ^tr {a 2 G^(e)a 2 G A (e - 
uj)), and three triplet (L = 1,2,3) channels, Pq M = 
^tr{a L a 2 G^(e)(T2(r M G A (e - <*>)), as g wl oc P£° - 
Em=i,2,3 P c /M [ Ref - 0]- In the absence of SO cou- 
pling and Zeeman splitting, Cooperons Pq M split into 
completely independent channels: one singlet and three 
triplet, P = 5P, where § = S LM and P obeys the dif- 
fusion equation. The SO coupling and Zeeman splitting 
mix up various components B and split their spectra, 
which modifies the diffusion equation into the matrix 
equations TLP(X, X') = 5 ■ 5(X, X'), 

f[ = 7 <5 + ie z fj - D(Sd Xl + iS2A x - ^A" 1 ) 2 (15) 



-D{Sdx 2 +iS2A 2 + iS 1 X. 
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where S^ M ■ 

-KLM 



-IE 



KLM 



are spin-1 operators (K 
1,2,3), e KLM is the antisymmetric tensor (K,L,M 



1,2,3), and rj LM = IlSom + SqlIm indicating that co- 
herence between oppositely polarized electrons is lost on 
the time scale of e% , and D is the classical diffusion 
coefficient. Equation ( J15| ) is supplemented with the the 
boundary condition at the edge of the dot characterised 
by the normal direction nu = (ni, 112), 



5{V + i2A) - i ni S 2 \i l + inzSiXz 1 P = 0. (16) 



The correspondence between random matrix theory 
description of a disordered system and diagrams is usu- 
ally transparent in the zero-dimensional (OD) approxi- 
mation in the diffusion problem, E^ — > 00, when the 
lowest modes are taken in the coordinate-independent 
form and coupling to higher modes is treated as a 
perturbation. Here, the boundary condition in Eq. 
( Jig ) requires the use of rotation to a local spin- 
coordinate system, P = OPO^ 1 , prior to making 
the OD approximation, with O — exp{i[SxX2\ 2 — 



S , 2XiA]~ 1 ]}exp{— iip B (X)S3}exp{— iipA^X)} where har- 
monic functions <p transform the symmetric gauge in 
Eq. (|j) to such a gauge, where vector potentials on the 
boundary are tangential to it. This eliminates the lowest 
orders SO coupling terms from the boundary condition, 



and, in a small dot Q L^ 2 <C Ai,2, can be followed by a 
perturbative analysis of extra terms generated by rota- 
tion O in Eqs. dig). This step results in the OD matrix 
equation for the Cooperon, 
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The form of Eq. (fLTJ) is applicable beyond the diffu- 
sive approximation as it follows from purely the symme- 
try considerations. The difference in the third term in 
brackets reflects the addition or subtraction of the Berry 
and Aharonov-Bohm phases, as was pointed out in Ref. 
fllPf . The expression for the weak localization correction 

can be found from Eq. ( |i"7J ) as g w \ ex tr {P[5 — S 2 ]}. In a 
dot with t b x = and 7, e Z) D/\\ 2 < E T 0, this yields: 
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where we used the fact that e\ < ez. It is interesting 
to notice that the application of the Zeeman field alone 
(in plane magnetic field Q) does not suppress the weak 
localization completely at ez <C Et . In the opposite case 
of ez > Et, it does (though has to be studied beyond the 
universal limit). However, the effect of such a strong in- 
plane fields on orbital motion becomes already sufficient 
to suppress the weak localization |ll[ . 

The form of Eq. (Q) in the experimentally easier 
achievable 'high'-energy crossover 1^4^7 is limited by 
only two first terms and suggests a possible procedure 
for measuring the ratio A1/A2. By fitting experimen- 
tal magnetoresistance data to g w i(ez) in Eq. (|l8|), one 
would determine the characteristic in-plane field B. For 
a dot with a strongly anisotropic shape, such a parameter 
would depend on the orientation of an in-plane magnetic 
field. In particular, B can be measured for two orienta- 
tions of B = Bl: namely, Smo] f° r l = [HO] and £>[iio] for 
I = [110]. One should also make a simultaneous measure- 
ment of two characteristic fields £>f 110 ] ano - ^filoi ^ n a ^°^ 
produced on the same chip by rotating the same litho- 
graphic mask by 90°. The anisotropy of the SO coupling 
can be then obtained directly from the ratio 



e [iio]S[ii ]/B[ii ]S[ lio] J = (A1/A2) 



independently of the details of sample geometry. 

The other interesting feature may be observed in the 
weak localisation g w \ in the regime of a crossover 2^2u 
driven by a weak perpendicular magnetic field for ez = 0. 
Since the SO coupling effect in Eqs. (§|,[17|) acts as a 
homogeneous magnetic field distinguishing between up- 
and down-spin electrons, the external field can be used 



to compensate the effect of the SO coupling for one spin 
component ]10| , which would produce a dip in the weak 
localisation correction. Indeed, at 7 C e™, the function 
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(7 + r B 1 )(7 + r B 1 +2ef) 
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has a minimum at the value of the field 2eB z — 
cfi./(AiA 2 ), independently of the sample geometry, which 
should provide a very accurate measurement of A1A2. 
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